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ABSTRACT 

We present a R ~ 10,000 M-band spectrum of LLN 19 (IRAS 08470-4321), a heavily em- 
bedded intermediate-mass young stellar object located in the Vela Molecular Cloud, obtained 
with VLT-ISAAC. The data were fitted by a 2-slab cold-hot model and a wind model. The 
spectrum exhibits deep broad ro-vibrational absorption lines of 12 CO v = 1 <— and 13 CO 
v = 1 4— 0. A weak CO ice feature at 4.67 /im is also detected. Differences in velocity 
indicate that the warm gas is distinct from the cold millimeter emitting gas, which may be 
associated with the absorption by cooler gas (45 K). The outflowing warm gas at 300^1-00 K 
and with a mass-loss rate varying between 0.48 x 10~ 7 and 4.2 x 10~ 7 M© yr _1 can ex- 
plain most of the absorption. Several absorption lines were spectrally resolved in subsequent 
spectra obtained with the VLT-CRIRES instrument. Multiple absorption substructures in the 
high-resolution (R= 100,000) spectra indicate that the mass-loss is episodic with at least two 
major events that occurred recently (< 28 years). The discrete mass-loss events together with 
the large turbulent width of the gas (dw=10-12 km s _1 ) are consistent with the predictions of 
the Jet-Bow shock outflow and the wide-angle wind model. The CO gas/solid column density 
ratio of 20-100 in the line-of-sight confirms that the circumstellar environment of LLN 19 is 
warm. We also derive a 12 C/ 13 C ratio of 67 ± 3, consistent with previous measurements in 
local molecular clouds but not with the higher ratios found in the envelope of other young 
stellar objects. 
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1 INTRODUCTION 

The gas dynamics of the inner surroundings of low- and high-mass 
young stellar objects (YSOs) is a complex interplay between ac- 
cretio n, settling and ejection (e.g. JStahler & Palla 200^; lArce et al.l 
120071) . The gas and dust of the collapsing envelope are accreting 
onto the YSOs or settling onto discs in Keplerian rotation. Mean- 
while, powerful jets and winds remove the excess angular momen- 
tum dispersing the envelope and allowing the system to contract. 
The release of gravitational energy and the radiation of the YSOs 
heat the inner circumstellar gas to a few hundred to a few thousand 
K. Mid-infrared observations best probe these regions for two rea- 
sons. First, dust grains absorb and scatter less in the infrared than in 
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the visible. Second, atoms and molecules at a few hundred Kelvin 
emit preferably in the mid-infrared. 

Molecular ro-vibrational transitions trace the large range in 
density, temperature and ultraviolet radiation field encountered in 
circumstellar discs and envelopes. While rotational lines in the mil- 
limeter wavelength range trace the cool outer disc and envelope 
(T < 100 K), the mid-infrared fundamental emission (v = 1 — >■ 0) 
and absorption (v — 1 4— 0) probe the warm part (T= 100-900 K) 
and the overtone emissions ( Av=2) in the near-infrared the ho t dust 
free region at T >1500 K dBik & Thill2004l : iThi et alj|2005h . CO 
ro-vibrational lines seen in absorption against the continuum gen- 
erated by warm dust can probe cold as well as warm gases along 
the line-of-sight. 

Fundamental absorption lines of CO, the second most abun- 
dant molecul e after Fb, are detected in high-mass yo ung stellar 
objects (e.g., iMitchell etai]|l99Cl : IScoville et al1 fl983). The line 
profiles indicate the presence of quiescent cold gas located in the 
envelope as well as outflowing warm gases. Ro-vibrational CO 
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gas lines in emissi on and in absorptio n have also been seen in 



the envelope (e.g. iBoogert et alj 120021 : IPontoppidan et al 
ll l2007l) " 



2003; Bri ttain et al J 2007) or in the inner disc iLahuis et al 



200: 



200l 



of low-mass embedded young stellar objects. Be sides CO, ab- 
sorption lin es of formaldehyd e dRoueff et al.ll2006h . C 2 H 2 , HCN, 
PCS. NHa jEvans et aljl99lh . "C 12 CH 2 , CH 3 , NH 3 , HNCO, CS 
dKnez et al. 2009), and methane dBoogert et alj2004h are observed 
towards high-mass YSOs with high-dispersion spectrometers from 
the ground. 

In this paper, we present an analysis of the CO funda- 
mental absorption lines toward LLN 19 (IRAS 08470-4321), an 
intermediate-mass YSO located in the Vela Molecula r cloud D at 
an estimated distance of 700 pc dLiseau et alj|l992h . The lumi- 
nosity of the object (L t = 1600 Lq) corresponds to a ~ 6-7 Mq 
protostar JPalla&Stahlerlll993h . LLN 19 is a deeply embedded 
(Ay ~ 45 mag) Y SO with moderat e water and CO ice abundance 
in its line-of-sight JThi et alj|2006h . Millimeter continuum emis- 
sion shows that LLN 19 is surrounded by 3.5 Mq of gas and dust 
in a 24 " beam. Part of t his matter likely fo rms a disc and the rest 
remains in the envelope dMassi et alj |"l999). A molecular outflow 
is detected in the 12 CO J = 1 ->■ 0, C 17 Q J = 1 -> and CS 
J = 2 — > 1 maps of LLN 19 dWouterloot & Brandl 19991) . The total 
mass of the out flow reaches 55 Mq if a kinematic distance of 2.24 
kpc is adopted dWouterloot & BranJl9"99^ . lLorenzetti et al.1 fc002h 
found shocked H 2 emission around LLN 19 but could not defini- 
tively show that LLN 19 is the source of the outflow. Observation 
evidence for the episodic nature of large scale molecular outflows 
is mo unting in low- and high-mass YSOs (e.g., lArce& Goodman 
l200lh . However, evidence for episodic warm inner winds exists 
primarily for high-m ass YSOs from CO absorption studies (e.g., 
Mitc hell et all [199 1). Our work extends the study of inner warm 
winds to lower luminosity objects. 

We describe the observation and data reduction procedures in 
Sect. [2] and give a first analysis of the results in Sect. [3] We ana- 
lyze the data using a wind model for the ISAAC data and a wind- 
envelope model for the CRIRES data in Sect. [5] and Sect. [6] and 
discuss the implications of our analysis on the gas and ice content 
in the environment of LLN 19 in Sect. [7] Finally, we conclude in 
Sect.® 



flat-fielding, distortion correction, etc.). Because of the scarcity of 
halogen lamp lines in the Af-band, the spectrum was wavelength 
calibrated by comparing the strong atmospheric absorption lines 
with high signal-to-noise ratio spectra of the atmosphere above 
Paranal measured by the ESO staff. The uncertainties in transmis- 
sion are do minated by systematics from the te lluric line removal 
procedure jPontoppidan & v an Dishoeck 2005). The precision of 
the wavelength calibration is ~ 5 km s" 1 , i.e. 1/6 of the resolution 
dv=30 km s _1 , over the entire sp ectrum. Details on t he dat a reduc- 
tion procedures are described in Pontoppid an et al.l d2003l) . At the 
time of the observation (7/1/2001), u obs =VLSR(LLN 19)+4.2 km 

Subsequent high-resolution observations at R= 100,000 (dv=3 
km s _1 ) of LLN 19 were obtained with the Cryogenic Infrared 
Echelle Spectrograph (CRIRES) at the VLT on April 26 th 2007. 
CRIRES provides a resolving power of up to 10 in the spectral 
range from 1 to 5[im when used with a 0.2 "slit. The simultaneous 
spectral coverage between 4.64 and 4.85 /im was achieved using 
a mosaic of four Aladdin III InSb arrays, which combine into an 
effective 4096 x 512 focal plane detector array in the focal plane. 
We used the Adaptive Optics system MACAO (Multi- Applications 
Curvature Adaptive optics) to optimize the signal-to-noise ratio and 
the spatial resolution. In our observation, the signal-to-noise ra- 
tio reaches ~ 180. As for the ISAAC observations, chopping and 
nodding along the slit were used and standard stars were observed 
at the same airmass to correct for telluric absorption lines. The 
CRIRES data were reduced using standard procedures that includes 
flat-field correction, detector non-linearity correction, and lineari- 
sation of the spectra in both the dispersion- and cross-dispersion 
directions. We used a telluric model to wavelength-calibrate the 
spectrum and to blank-out the telluric absorptions after division by 
the standard. The width of the continuum is 0.35 ", consistent with 
the AO-corrected spatial resolution. The CRIRES data are part of 
a larg e programme that observ ed ~100 YSOs and protoplanetary 
discs (Pontoppidan et al. 2008). A handful of the observed YSOs 
show high velocity CO absorption features that can be ascribed to 
winds. In particular, the spectrum of LLN 19 has the deepest CO 
absorption features. 



2 OBSERVATION AND DATA REDUCTION 

We observed LLN 19 (IRAS 08470-4321) at a spectral resolving 
power of R w 10, 000 (resolution of 30 km s" 1 ) by setting the slit 
width to 0.3"with the Infrared Spectrometer And Array Camera 
(ISAAC) mounted on the Very Large Telescope-Mil (VLT-UT1) 
as one of the targets of a ESO large p rogramme to study the ic e 
and gas around young stellar objects ( Ivan Dishoeck et aT]|2003b . 
The slit width corresponds to a physical distance of 210 AU at the 
distance o f LLN 19 (700 pc). We adopted the distance of 700 pc 
derived bv lLiseauet"aIUl992h over the kinematic distance of 2.24 
kpc. The closer distance is consistent with the YSO LLN 19 be- 
ing a member of the Vela Molecular Cloud. We obtained a M-band 
spectrum from 4.5 to 4.8 \iva of LLN 19. We also observed standard 
stars at same airmass with the same instrument settings. We restrict 
our analysis to the 4.55-4.8 /im region because the signal-to-noise 
ratio is too low below 4.55 /im due to strong atmospheric C0 2 ab- 
sorption. The detector for 3-5 [im observations is a 1024 x 1024 
Aladdin InSb array. We reduced the data using an in-house software 
written in the Interactive Data Language (IDE). The data reduction 
procedure is standard for infrared observations (bad pixel maps, 



3 OBSERVATIONAL RESULTS 

The continuum subtracted ISAAC spectrum in Fig.Q]is dominated 
by deep gas-phase 12 CO ro-vibrational v = 1 •<— P-branch (J + 
1 J' = J+ 1, J" = J) and P-branch (J + 1 — > J) absorption 
lines. A shallow absorption feature due to solid CO is also present 
in the spectrum. 

12 CO R(0) and R(l), and P-branch transitions from 1 to 13 
have been observed with CRIRES and a sample of those absorption 
lines is shown in Fig. [2] The absorption lines show multiple sub- 
structures, which correspond to the blue wing absorption seen in 
the lower-resolution ISAAC spectrum (Fig. [3J. Such data are avail- 
able only for the lower J lines, probing the cooler gas. 

3.1 CO ice feature 

A weak feature centred around 4.675 fim (2140 cm -1 ) correspond- 
ing to CO-ice is detected. The total solid CO column density is (4.3 
± 0.6) x 10 16 cnT 2 using the band strength of A(CO)= 1.1 x 
10~ 17 cmmolec" 1 at 14 K jjiang et al.ll 19751 : ISchmitt et alJI 19891: 
iGerakines et al.l ll995) to convert the integrated absorption to col- 
umn density. The detailed analysis of the CO ice components in 
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Figure 1. VLT-ISAAC M-band spectrum of LLN 19. The positions of the 12 CO and 13 CO lines are labelled. Wavelength regions with poor atmospheric 
transmission are omitted and labelled with a ffi symbol. 



relatio n to other YSOs in the Vela cloud is provided by iThi et"al1 
J2006h . 

3.2 Gaseous components 

3.2.1 ISAAC data 

We first focus on the ISAAC spectrum. The LLN 19 spectrum is 
particularly rich in gas-phase lines. We detect the hydrogen recom- 
bination emission line Pi/3, which is an indicator of accretion/wind 
activity and dense gases (n c ~ 10 8 cm -3 ), at 4.653 /im. The broad 
wings (AV=300-500 km s _1 ) of the line can be ascribed to a stel- 
lar/disc wind, especial ly in low- and interm ediate-mass stars sur- 
rounded by discs (e.g. ICarmona et al.ll2005h . For high-mass stars, 
the stellar wind is driven by radiative pressure exerted on the ionic 
lines. We tentatively detect the pure-rotational H2 — 5(9) line 
at 4.695 fim with an equivalent width of (27 ± 12) x 10" 3 cm -1 
(la). The line is also seen in the CRIRES data (see the lower-right 
panel of Fig. [2} with an equivalent width of (11.4 ± 3.5) x 10 -3 
cm -1 . The detection of the H2 line in emission, whose energy level 
lies at 5005.73 K, testifies to the presence of hot molecular gas. 

12 CO R— and P-branch lines, up to J=10 and J=19 respec- 
tively, are clearly seen (Fig. [TJ. The lines are up to 60 % deep at 
the spectral resolution of 30 km s _1 . Columns 1 and 2 of Table Q] 



contain the rest laboratory and observed wavenumbers. The quies- 
cent cold gas around LLN 19 is located at VLSR(mm) = 12 km s - 
from millimeter observations of 12 CO J = 1 — > in a 46" 
beam dWouterloot & Brand fl99^) . By contrast, the CO gas seen 
in the mid-infrared is moving by an average ulsr = -25 km s _1 
(± 5 km s _1 ). The line velocities are not randomly scattered around 
the average velocity shift value. The low-J lines have line widths 
smaller than the high-J lines, suggesting that quiescent cold gas 
may also be present in the line-of-sight, but the spectral resolution 
is too low to resolve it from the outflowing gas. The line profiles are 
not symmetric and can be decomposed into two Gaussian compo- 
nents: a deep component centred at ulsr = -25 =F 5 km s -1 and Full 
Width at Half Maximum FWHM= 33 km s" 1 and a broad shal- 
lower component at «lsr = -75 =p 5 km s~ with FWHM = 59 
km s" 1 . Fig. [3] displays a normalised absorption line profile aver- 
aged over all the lines and its decomposition into the two compo- 
nents. The lines do not broaden with increasing quantum rotational 
number for the transition. The decomposition is arbitrary but allows 
a simple analysis of the two components. 

In the rest of the paper, we will refer to the main (at -25 
km s -1 ) and blue (at -75 km s -1 ) component. At the resolu- 
tion of ~ 30 km s" 1 , the main and blue CO absorption lines are 
marginally resolved; the depth of the absorption lines indicates 
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Figure 2. Absorption lines in open circles towards LLN 19 observed with VLT-CRIRES at a resolution of 3 km s~ 1 . 
transmission are omitted. The long-dashed lines are the substructure fits and the continuous lines are the total fits, 
components. The lower-right panel shows the tentatively detected H2 S(9) line at 4.6947 fim by VLT-CRIRES. 



4.6940 4.6945 4.6950 4.6955 
Wavelength (micron) 

Wavelength regions with poor atmospheric 
which are the sum of all the Gaussian sub- 



that the intrinsic width Ai; is ~ VFWHM 2 - Resolution 2 = 
V33 2 - 30 2 = 13.7 kms" 1 if the main absorption profile 
is caused by a single component. High velocity blueshifted gas 
is commonly found to wards high-mass dMitchelletalJll99lh and 
low-mass YSOs (e.g jLahuis et al]2006h . 



Weaker CO i?-branch absorption lines are detected up to 
J =22 as well and are labelled in Fig.Q] The blue part of the average 
15 CO profile is contaminated by other features, mostly by the wings 
of the 12 CO feature. The 13 CO lines are seen at the same velocity 
than their 12 CO counterpart. However, the weakness of the 13 CO 
lines forbids an accurate determination of the shifts. 



The equivalent widths of the main and blue component are 
given in Tables Q] and [2] The normalised equivalent width ratios 
Wo/v = W\/X between 12 CO and 13 CO lines increase with in- 
creasing rotational quantum number but remain much lower than 
the interstellar isotopic value of 12 C/ 13 C = 60-65, suggesting that 
the 12 CO lines are optically thick. 



3.2.2 CRIRESdata 

The CRIRES spectra show that the two components seen in the 
ISAAC spectrum actually consist of the absorption by four clumps 
at 9, -1 1.5, -27, -43.5, and -70 kms" 1 . The components at 9, -1 1.5, 
and -27 km s _1 correspond to the main component in the ISAAC 
data while the -43.5 and -70 km s _1 components correspond to the 
blue wing. The absorption at 9 km s~ 1 may be saturated at the reso- 
lution of 3 km s -1 . As a result the measured widths of the 9 km s _1 
feature may be broader than the intrinsic widths. 

We modelled the absorption lines by a combination of four 
Gaussians with fixed velocities Vi (9, -11.5, -27, -43.5, and - 
70 km s _1 ). We performed our Gaussian decomposition in the 
transmission domain since we are interested in the equivalent 
widths of each component. The fitting function has the form: 

i 

where Ai and 6vi, the widths and amplitudes of the Gaussians, are 
found by the fitting routine. The width and amplitudes of the Gaus- 
sians are free parameters of the fits. We used the Shuffle Complex 
Evolution optimization code (see Sect. I5.lt to determine the best 
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Figure 4. Gaussian width (left panel) and equivalent width (right panel) for each individual substructure and each transitions (x axis) for the CRIRES. The 
number in each panel indicates the velocity of the substructures. The main causes of uncertainties are missing data due to strong telluric absorptions. 
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Figure 3. Comparison of the 12 CO and 13 CO average ISAAC profiles. The 
profiles are averages of all the detected lines to increase the signal-to-noise 
ratio. The profiles are corrected for the Earth velocity at the day of the obser- 
vation and are normalised to an absorption depth of 0.4 to ease the compar- 
ison. The model decompositions are shifted by +0.2 for clarity. Both mean 
profiles were shifted by -9 km s — 1 . 



fits. The fit for each individual line is shown in Fig. [2] whereas the 
FWHM of the Gaussians and the equivalent widths are shown in 
Fig. EH 

The largest equivalent widths {EW > 0.1 cm -1 ) correspond 
to the outflowing component at -1 1 km s~ , The quiescent envelope 
component at 9 km s~ shows decreasing equivalent widths with 
increasing J values, testifying of a cold component. The large EW 
for P(13) in the CRIRES data for the component at 9 km s _1 is 



most likely an artifact caused by missing data. The gas equivalent 
widths do not show similar clear trend with increasing J values. 

The average FWHM of the Gaussians is consistent with an 
intrinsic turbulent width of ~ 13.7 km s _1 inferred from the ISAAC 
data. The parameters will be analysed in the following sections. 



3.2.3 Comparison between ISAAC and CRIRES data 

The measured equivalent widths derived from ISAAC and CRIRES 
data are compared in Fig. [5] Discrepancies are small for most tran- 
sitions (<15%) except for P(l) and R(0). Generally the CRIRES 
equivalent widths are larger than the ISAAC equivalent widths. 
The lower resolution ISAAC data may suffer from systematic er- 
rors in the determination of the continuum, especially when the 
12 CO blue-wing absorptions reach the wavelengths of the 13 CO ab- 
sorption. The large differences for P(l) and R(0) may stem from 
the incomplete profiles in the CRIRES data-set at those transitions 
due to strong telluric absorptions. In case of strong telluric ab- 
sorptions, the data are interpolated to compensate for the miss- 
ing data for the determination of the equivalent width and this 
procedure can introduce significant systematic errors, in our case 
an overestimate of the actual absorption (see Fig. [2}- Strong tel- 
luric absorption features also introduce errors in the ISAAC equiv- 
alen t widths since the raw ISAAC data are ratioed by a standard 
star Jpontoppidan & van Dishoeck 2005). Alternatively the choices 
made in the fitting of the CO ice absorption feature affect the con- 
tinuum flux around the P(l) and R(0) lines more strongly in the 
ISAAC spectrum than in the CRIRES spectrum (see Fig.[TJ. There- 
fore the P(l) and R(0) equivalent widths measured by ISAAC may 
be underestimated. In summary, the differences in EW between the 
two datasets can probably be ascribed to systematic errors in inter- 
polating missing data and in cancelling telluric absorptions. 
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Figure 5. Plot of percentage difference versus mean value between equiv- 
alent widths derived from the ISAAC and the CRIRES observations. The 
systematic differences are below 20% except for the P(l) and P(l) transi- 
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Table 1. 12 CO v = 1 <— R- and P-branch line positions and equiva- 
lent widths of the middle component (mc) and blue component (be) toward 
LLN 19. 



Line 






W (mc) 






(cm 




(10~ 3 


cm -1 ) 


R(0) 


2147.081 


2147.19 


208.3 ± 1.6 


61.4 ±0.3 


R(l) 


2150.856 


2150.96 


263.4 ± 14.3 


77.7 ± 3.0 


R(2) 


2154.596 


2154.72 


279.7 ± 10.9 


82.3 ± 2.3 


R(3) 


2158.300 


2158.43 


300.5 ± 2.5 


88.7 ± 0.5 


P(4) 


2161.969 


2162.07 


306.9 ± 0.5 


139.9 ± 0.2 


R(5) 


2165.601 


2165.74 


296.8 ± 0.7 


135.2 ± 0.2 


R(6) 


2169.198 


2169.34 


292.5 ± 6.9 


133.3 ±2.3 


R(7) 


2172.759 


2172.87 


273.5 ± 5.2 


145.9 ± 2.0 


R(8) 


2176.284 


2176.42 


281.3 ±7.3 


171.4 ±3.2 


R(9) 


2179.772 


2179.93 


279.4 ± 8.9 


123.6 ± 2.8 


P(10) 


2183.274 


2183.38 


275.6 ± 4.8 


81.1 ± 1.0 


fl(ll) 


2186.639 


2186.84 


291.9 ± 6.3 


85.9 ± 1.3 


R(12) 


2190.010 


2190.24 


271.9 ± 1.2 


80.0 ± 0.3 


R(13) 


2193.359 


2193.66 


192.0 ± 1.2 


70.6 ± 0.3 


R(14) 


2196.664 


2196.88 


239.8 ± 1.0 


70.6 ± 0.2 


P(15) 


2199.931 


2200.03 


232.4 ± 3.9 


191.6 ±2.3 


P(16) 


2203.161 


2203.42 


231.6 ±5.3 


107.9 ± 1.8 


R(17) 


2206.354 


2206.53 


137.5 ±27.4 


149.5 ±21.3 


R(18) 


2209.509 


2209.80 


213.9 ±21.4 


133.2 ± 9.5 


P(19) 


2212.626 


2212.94 


218.9 ± 31.9 


96.6 ± 10.1 


P(l) 


2139.426 


2139.59 


182.0 ± 9.6 


53.7 ± 2.0 


P(2) 


2135.546 


2135.77 


255.3 ± 1.9 


75.3 ± 0.4 


P(3) 


2131.632 


2131.74 


272.3 ± 1.3 


87.5 ± 0.3 


P(4) 


2127.684 


2127.84 


277.8 ± 3.5 


102.4 ± 0.9 


P(5) 


2123.699 


2123.86 


278.4 ± 1.8 


108.6 ± 0.5 


P(6) 


2119.681 


2119.86 


286.9 ± 2.5 


135.4 ±0.8 


P(7) 


2115.629 


2115.78 


267.3 ±2.1 


138.0 ±0.8 


P(8) 


2111.543 


2111.71 


255.7 ± 12.0 


124.6 ± 4.2 


P(9) 


2107.424 


2107.57 


229.7 ± 6.7 


105.0 ± 2.2 


P(10) 


2103.270 


2103.40 


238.1 ± 3.7 


108.8 ± 1.2 



3 <r fitting error. 



Table 2. 13 CO v 
toward LLN 19 



1 <— P-branch line positions and equivalent widths 



Line 




v obs (mc) 


W^rnc) 


Wp(bc) 






) 


(10 — 3 cm 


/ 


P(2) 


2106.903 


2106.82 


42.5 ± 4.5 a 


<5.9 


P(3) 


2110.447 


2110.42 


48.5 ± 4.5 




P(4) 


2113.958 


2113.85 


39.0 ±4.5 


-9.6 


R(5) 


2117.436 


2117.57 


34.4 ± 4.5 




P(10) 


2134.318 


2134.38 


27.4 ± 4.5 


-21.7 


P(12) 


2140.833 


2141.05 


21.3 ±4.5 




P(13) 


2144.039 


2144.15 


19.9 ±4.5 




P(16) 


2153.447 


2153.59 


19.7 ±4.5 




P(17) 


2156.514 


2156.45 


22.4 ± 4.5 




P(18) 


2159.545 


2159.83 


17.7 ± 4.5 




P(22) 


2171.317 


2171.75 


16.4 ±4.5 





a formal 3 cr error dominated by the flux calibration accuracy. 

4 EQUIVALENT WIDTH ANALYSIS OF THE MAIN 
COMPONENT USING ISAAC DATA 

4.1 Curve-of-growth analysis 

We performed a classical curve of growth analysis of the 
ISAAC data, which relates the measured equivalent width Wt, (in 
cm -1 ) at wavenumber v (in cm -1 ) to the CO column density Njn 
(in cm" ) in the initial lower rotational state J" assuming a Gaus- 
sian profile with an intrinsic line width called the Doppler parame- 
ter bo = Avf W hm/2 log(2) 1//2 , w here A« F whm is the Full With 
at Half Maximum in km s" 1 (e.g.. lSpitzei]ll978l . Eq. 3.22), and a 
slab geometry. We will focus on the ISAAC main component. The 
equivalent widths Wp and Njn are linearly proportional only in 
the optically thin regime: 



s 



.7" 



fj 



■„Nj 



(2) 



where fji <_,/// is the oscillator streng th (dimensionless) of the 
transition jGoorvitch~& Chacker iaiil[l994f) . I„ is the background 
intensity and Sjn is the transition source function in the extended 
region covered by the beam and the other symbols have their usual 
meaning. We assume here that Sjn /I u <SC 1. Therefore the derived 
column densities are lower limits to the actual values. 

Transitions that were detected in both 12 CO and 13 CO 
with good signal-to-noise ratio in the ISAAC data are shown in 
Fig. [6] The vertical axis corresponds to the observed Wp/v, 
while the horizontal axis corresponds to the different values of 
Nj„( 13 CO)fj,^j„( 13 CO) lu and Nj,,{ X2 CO)fj,^j„(™CO)lv. 
Theoretical curve-of-growths computed assuming various values 
for the Doppler parameter bo are overplotted in dashed lines. The 
column density in a given level J" is computed by the equation 

Nj,,( 13 CO) = b D r( 13 CO)i>/(1.497 x 10" 2 fj.^jn). (3) 

Once the r( 13 CO) values are derived from the relationship be- 
tween t and Wc for each v alue of bp via the theoretical function 
F(t) = Q.^c{Wc/v)/bu jSpitzer|[l978l) . the column density of 
12 CO at level J" is 



Njn ( 12 CO) = <( 12 C/ 13 C)/Vj»( 3 CO). 



(4) 



The ratio between the isotopologues 12 CO/ 13 CO, C( 12 C/ 13 C), is 
an adjustable parameter. In Fig. [6] the 13 CO data lie on the linear 
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Figure 6. Curve of growth. Only the 12 CO equivalent widths (EW) are af- 
fected by saturation and thus can be used to constrain the Doppler parameter 
6 D . The 13 CO EW are on the linear part of the curve. The fit for the 13 CO 
EW are also shown in the inset for better clarity. The best fit curve to the 
data is obtained with b n =ll ± 1.5 km s _1 and 12 C/ 13 C = 67 ± 3. 



part of the curves; therefore, the lines are optically thin whereas the 
12 CO lines lie on the saturated part of the curve and are therefore 
optically thick. 

The best fit parameters were found by minimizing the differ- 
ence between the data and the empirical curve (least squares). First 
the £( 12 C/ 13 C) parameter was varied. Varying £ will only move the 
12 CO points left or right in the curve-of-growth diagram. For each £ 
value x 2 values between the observed EWs and the synthetic curves 
were determined by varying br>. The determination of £( 12 C/ 13 C) 
and 6d is independent on the actual CO column density and ex- 
citation. The best fit were obtained for <~( 12 C/ 13 C)=67 ± 3 (1 <j) 
and 6 D =1 1 ± 1.5 km s" 1 (1 a). Since the EW derived from ISAAC 
are strictly lower limits, we should consider that £ = 67 is also a 
lower limit. The errors were determined from the \ 2 surface. The 
C( 12 C/ 13 C) value derived here is an average between the £ value in 
the warm wind and that in the cold foreground cloud. The CRIRES 
spectra show that the main component profile is a combination of 
multiple Gaussians. The £ value of each individual velocity com- 
ponents can be derived using CRIRES data in future studies. The 
average value of C( 12 C/ 13 C) of 67 ± 3 is close to the local inter- 
stellar value of 65-77 (e.g jBensch et aLkOOllTsheffer et al.l2007l : 
IWilson&Ro"odlll994l) . 



4.2 Rotational diagram 

4.2.1 Main component 

We used the optically thin 13 CO data to construct the rotational 
diagram based on the ISAAC data that is shown in Fig. UJ for the 
main component. In a rotational diagram for i?-branch absorp- 
tion lines the natural logarithm of the lower level column density 
of a transition J' <— J" divided by the degeneracy (2J" + 1), 




400 600 

E,../k (K) 



1000 



Figure 7. Rotational diagram of 13 CO for the main component. The fit by 
the 2-slab model is overplotted. 



(log [Njn / (2 J" + 1)]) is plotted against the lower level energy 
Ejtt /k (in K). If the level population is thermalised (LTE), the col- 
umn density at a given level J" is expressed as: 

exp(-Ej,,/Tk in ) 



Njn = iV(CO)(2J" + 1) 



(5) 



where iV(CO) is the column density of CO, Ejn is the energy of 
the J"th rotational level expressed in K, Q(Tkm) — kThinj Bq + 
1/3 is a good approximation of the partition function at the gas 
temperature T)y n (equals to the gas excitation temperature at high 
temperatures and densities) and Bo is the rotational constant with 
B /k= 2.765 and 2.644 K for 12 CO and 13 CO respectively. The 
data points in a rotational diagram can be fitted by a straight line 
for an isothermal medium at LTE. The gas temperature is then 
given by the negative of the inverse of the slope of the line. How- 
ever, fits with a single temperature gas to the observed values pro- 
vide bad results (Fig. [7}. Possible non- LTE population was c hecked 
using the simple analytical results of Mc Kee et alj dl982l) . High- 
J levels in a gas at non-LTE are underpopulated compared to a 
gas at LTE, which is not seen in our data. Instead, the data are 
best fitted with the combination of a cold and hot gas at LTE, 
similar to most young stellar objects where CO fundamental ab- 
sorption lines have been observed. The low- J absorption lines are 
formed in a cold region (T gas = 45 K, iV( 13 CO)~ 3.3 x 10 16 cm" 2 , 
iV( 12 CO)= 67 x iV( 13 CO) ~ 2.2 x 10 18 cm -2 ) while the high- J 
lines (J ^10) are formed in warm gas (T gas = 753 K, iV( 13 CO)~ 
10 17 cm -2 , iV( 12 CO)~ 6.7 x 10 18 cm -2 ). The results are sum- 
marised in Table|3] The LTE level population and the non detection 
of high-velocity gas in emission in the millimeter indicate that the 
warm component arises in a low surface area, high density region 
(n,H >10 7 cm -3 ) with a very small filling factor in the millimeter. 



4.2.2 High velocity gas 

The blue absorption component in the ISAAC data was assumed 
optically thin, because the wing is absent in 13 CO, as testified by 
the 12 CO and 13 CO profiles in Fig. [3] We constructed a rotational 
diagram, displayed in Fig. [8] taking into account the transitions 
from levels detected in the R— branch and bzj=ll km s" 1 to cor- 
rect for optical thicknesses. Two gas components have been fitted 
to the data. Warm gas at ~1 16iSjg K and with column density of 

K and column 



2.0^1 x 10 1 ' cm" 2 and hot gas at ~ 1338 



-2298 
-518 
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Table 3. Column densities and temperatures derived from the rotational diagram analysis. The upper part of the table corresponds to the parameters derived 
using the ISAAC data while in the lower part of the table corresponds to the parameters driven using the CRIRES data. The formal errors are \a from the 
least squares fits. The 12 CO column densities of the cold and warm components are derived using the optically thin 13 CO data while the column densities in 
the blue wing were derived using the 12 CO data. The fits to the wing data are poor, which explains the large errors in the derived temperatures and column 
densities. The column densities derived from the CRIRES equivalent widths are corrected for optical depth effects. 



Component 


T 


AT( 12 CO) 


AT( 13 CO) 




(K) 


(cm" 2 ) 


(cm" 2 ) 






ISAAC data 




Cold main 


45 ±4 


(2.2 ± 0.6) x 10 18 a 


(3.3 ±0.1) x 10 16 


Warm main 


753 ± 150 


(6.7 ± 0.6) x 10 18 a 


(1.0 ±0.1) x 10 17 


Cold blue wing 
Warm blue wing 


1338± 22 g 8 


2.0 +1 l x 1017 
3.3±5;| x 10 17 








CRIRES data 




at +9 km s _1 
at-11.5kms _1 
at -27 km s — 1 
at -43.5 km s -1 
at -70 km s -1 


54+? 
014+33 

OQ7+1050 
•"'-192 
9Q1 +66 

387± 7 f 


lA+^l x 10 17 
5.2+^' 8 x 10 17 
6.9+^| x 10 16 
4.7+^3 x 10 17 
1.9+.?; 7 x 10 17 





a from "CO using 12 CO/ ld CO=67. 



ln[N/(2J+1)] 



R-branch 
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Figure 8. Rotational diagram of the blue- wing component of CO. 

density of 3.3 1 \% x 10 17 cm -1 provide the best fit (Tabled. The 
total column density of high- velocity gas is a factor 10 smaller than 
that of the main absorbing gas using CO, assuming that the 12 CO 
lines are optically thin. 

4.2.3 CRIRES components 

We drew the rotational diagram for 1 CO for the individual sub- 
structures in Fig.[9]and the results are summarized in the lower part 
of Table[3]. The CRIRES lines probe the cool gas towards LLN 19 
with substructure temperatures ranging from 47 to 396 K. 

High-J 13 CO transitions up to R(22) were available in the 
ISAAC data set to constrain the high-temperature gas compo- 
nent but they were not observed with CRIRES. The lack of high- 
temperature gas tracers explains the lower overall gas temperature 



derived from the CRIRES data. Optically thick 12 CO lines were 
used to derive the column densities with the optical depth correc- 
tions that were applied to the blue- wing ISAAC data (bo = 11 
km s" 1 ). The total column density for the main component, which 
comprises the components at +9, -11.5, and -27 km s , is 7.3 x 
10 17 cm -2 , a factor ~10 lower than the value derived from the 
13 CO ISAAC data. The column density of the component at -27 
km s _1 is uncertain by a factor ~25. The column densities de- 
rived from CRIRES data have large systematic errors because of the 
assumptions made for the optical depth corrections (the Doppler 
width bo and the number of absorption features). In the ISAAC 
dataset the column density of the main component is dominated by 
the column density of warm gas. Unfortunately the warm gas can- 
not be probed in the CRIRES dataset, which is limited to transitions 
up to P(13) (see Fig. [2}. The systematic errors and the lack of high- 
J lines in the CRIRES dataset may explain the discrepancy in 12 CO 
column densities derived from the two datasets. The two velocity 
components at -43 and -70 km s _1 in the CRIRES data match the 
ISAAC blue-wing component well in terms of temperature (~ 290- 
390 K) and column density (1.9-4.7) x 10 17 cm -2 . The reason for 
the similarity may stem from the fact that the blue-wing parameters 
were both derived from optically thick 12 CO data. 



5 WIND MODEL 

In this section we focus on the model-fitting on the ISAAC data 
since they are more complete in the rotational level coverage 
than the CRIRES data. We used the velocity structure seen in the 
CRIRES data however, to bound the parameters. 



5.1 Model description and fitting 

The temperature and column density of the high- velocity gas sug- 
gest that the absorbing gas may be located in the mixing layer be- 
tween the outflow and the envelope. The simple analysis of the 
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Figure 9. Rotational diagram of 12 CO for each substructure components in the CRIRES data. Equivalent width for P(12) and P(13) were not used in the 
analysis of the 9 km s _1 substructure because the absorption features are compromised by strong telluric absorption (see Fig. [2). 



ISAAC data has assumed an arbitrary distinction between a main 
component and a blueshifted one. The CRIRES data have shown 
that the outflow rate of gas is not continuous but may instead pro- 
ceed by episodes of strong mass-losses. However, a simple static 
model like the one we adopt here cannot handle dynamical effects, 
although the model can still give an average value for the mass loss 
rate. 

We fitted the 12 CO and 13 CO synthetic transmission spectrum 
of a wind model to the transmission spectrum of LLN 19. The 
model focuses on the molecular part of the outflow only. 

The spherical wind compon ent of the model is a simplified 
version of the parametric model of lChandler et alj l ll995l) . The wind 
is divided into concentric shells located at distances R from the star. 
The wind is launched from Ri n and the size of the spherical wind 
is Rout- We assume an accelerated wind with velocity of the form: 
R \ q 



(V max -V Rin )(l-^ +V Rl . 



(6) 



where Vn in is the wind velocity at Ri n and V ma x is the asymp- 
totic wind velocity. The index q controls the variation of the wind 
velocity with R. The wind velocity does not vary in time in our 
model. 



Assuming a constant mass-loss rate m (in Mq yr _1 ) through- 
out the wind, the spherical wind velocity, mass-loss rate and col- 
umn density are related by the mass conservation equation 



n( 12 CO) = 



47rmH 2 -R 2 V / u ,(l — cosi 



n( 12 CO) 
n(H 2 ) 



(7) 



where m,H 2 is the gas molecular mass (m,H 2 = 2.2 a.m.u.), 
n( 12 CO) /n(H2) is the 12 CO abundance relative to H2, and <f> is the 
semi-opening angle of the wind (for a spherical wind, (j> = 7r/2). 
For a spherical wind, a shell of thickness AR has a column density 
AN( 12 CO)(R)=n( 12 CO)(R)xAR and <j> = tt/2. We assumed 
that the wind is fully molecular and that n( 12 CO)/n(H 2 )=10~ 4 . 

The gas temperature was not computed self-consistently, in- 
stead we adopted a power-law T rot = T in (R/ Ri„)~ v . The gas 
temperature power-law index can take positive and negative val- 
ues. If the wind cools radiatively and adiabatically by expansion, 
then p is positive. The wind can also be heated by shocks and tur- 
bulence. In that case the value for p can be negative. The 12 CO and 
13 CO level population are at LTE. 

Spectra were computed by a simple line+continuum ray- 
tracing with T gaa = Tdust = Trot = T vib . The ray-tracing in- 
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tegrates the radiative transfer equ ation from shell i—1 to shell i for 
one ray t aking into account du st jDraine & Leell 19841) and CO gas 
opacities l lChandra etal.lll996h : 

= ~ V)e' Tv{%) + B v {T{i)) (1 - e~ T ^) . (8) 

We removed a low-order polynomial as continuum to the in- 
tensity spectra to obtain synthetic transmission spectra T moc j. 

The number of free parameters amounts to 8: 9=[Ri„, Rout, 
Vn in , Vmax, q, in, Ti n , p]- The intrinsic turbulent line width bo 
is set to 10 km s _1 and the 12 CO/ 13 CO ratio is assumed to be 65. 
The number of data points n is 1020 (1024 minus 2 points at each 
edge). The number of degrees of freedom (dof) is thus 1012. The 
high number of parameters means that degeneracy between param- 
eters is unavoidable. We also imposed physically-motivated lim- 
its to the parameters. Assuming an outer radius of 100-500 AU, 
we can estimate the mass-loss rate using equation [7] a CO column 
density of 6.7 x 10 18 cm -2 and a wind velocity of 50 km s _1 . We 
found that the mass-loss rate is of the order of 10 -7 Mq yr _1 . The 
mass-loss rate ranges from 10 -8 to 10 -6 Mq yr _1 during the fit- 
ting procedure. The wind acceleration parameter q was allowed to 
vary between and 3. The gas temperature can decrease (positive 
p values) or increase (negative p values) as the wind expands. 

We used a non-linear global optimisation code based on the 
Shuffled Complex Evolution (SCE) algorithm, which mixes evo- 
lutionary and simplex algorithm jDuan et al]|l993f) , to minimise 
the Residual Sum Squared (RSS) between observational data and 
model outputs 

F(6) = RSS = Y, (To^i) ~ T mod (v 3 ,9)f = ^ ff(6), (9) 

3=1 3=1 

where T ^ a (i>j) is the observed transmission depth at wavenum- 
ber Vj and T mo d(^j, 9) is the corresponding modeled transmission 
with the k = 11 parameters. F(9) is the residual sum squared 
(RSS) or merit function and is related to the chi-squared. If the 
noise a is uniform, by definition xS=X 2 ^of = (RSS/cr 2 )/dof, where 
dof is the degree of freedom. Minimizing the RSS is equivalent to 
minimizing \ 2 ■ We obtained a reduced xt of 0.4, which suggests 
that the noise level is overestimated in average (i.e., the S/N of 99 
is underestimated). 

The SCE is a hybrid algorithm that uses a Genetic algorithm 
with local p opulation evolution r eplaced by a Nelder-Mead Sim- 
plex search (Nelder&Mead 1965b . The Genetic algorithm part en- 
sures that a global minimum can be found while the Nelder-Mead 
simplex method provides fast local minimum searches. The SCE 
method is efficient in finding the global minima of non-linear prob- 
lems with multi local minima. 

5.2 Model results and error budget 

The best synthetic spectrum is compared to the observations in 
Fig. [10] After subtraction of the observed transmission by the wind 
model, the residual still show weak absorption in the J < 3 lines, 
testifying of the presence of cold (T < 100 K) material, maybe un- 
related to LLN 19 where water and CO ice are present. This cold 
component was analysed in Sect. 14.11 The fit also fails to explain 
the relatively noisy high- J absorption lines above v= 2200 cm -1 

For nonlinear models there is no exact way to obtain the co- 
variance matrix for the parameters. Instead we use the asymptotic 
approximation and provide the 95% confidence interval (see Ap- 
pendix for more details). The parameters are relatively well con- 
strained (see Table [4}. The asymptotic error method gives unreal- 



Table 4. Best fit parameters to the ISAAC data. No velocity correction has 
been applied. C. I. means confidence interval computed from the asymptotic 
method. The x 2 + 1 Confidence interval is also provided. 



Units Best fit 95% C. I. best < (x 2 + 1) C. I. 

value 





AU 


1.0 


±0.5 


0.1 


3.9 


Rout 


AU 


480 


±99 


74 


3668 


T 


K 


308 


± 10 


69 


925 


p 




-0.15 


±0.26 


-0.32 


0.40 


V Rin 


km s — 1 


16.1 


± 0.03 


14 


19 


Vmaa; 


km s — 1 


87 


±3 


11 


270 


q 




1.1 


±0.3 


0.01 


2.0 


ra 


l(r 7 Mq yr" 1 


4.2 


±0.7 


0.3 


11.6 



S/N 99 
Xl 0.4 



istically small errors on the parameters. In addition we provide the 
possible range for the parameters where xt < best xt + !• 

The best fit parameters for the wind model are summarised in 
Table [4] The wind base is at 1 AU and reaches 480 AU. Infrared 
absorption measurements do not constrain the spatial extent of the 
gas along the line-of-sight. Assuming gas densities greater than 10 6 
cm" 3 for the CO rotational level population to be at LTE and col- 
umn density of 5 x 10 22 cm -2 , the outflow radius should be smaller 
than 4500 AU, consistent with our derived value. The total 12 CO 
column density in the wind is 3.5 x 10 18 cm -2 , lower but compa- 
rable to the total of 9.9 x 10 18 cm -2 found using the simple two 
slab analysis. The wind temperature is 308 K at the base at 1 AU 
and it stays relatively constant through the warm wind with a slight 
increase (p=-0. 15). However, the temperature variation is also con- 
sistent with a temperature decrease within the errors. 

The asymptotic maximum wind velocity V ma x is 87 km s~ 
starting at 16 km s" 1 . The wind velocity index q is 1.1. The mass- 
loss rate is 4.2 x 10~ 7 M Q yr -1 . 



6 WIND-ENVELOPE MODEL 

In this section, we attempt to fit the CRIRES data with the same 
code but with added parameters because the low-J lines indicate 
that extra material at low velocity is present in the line-of-sight. In 
addition to the wind-model, we modelled simultaneously the ab- 
sorption by a quiescent slab component in addition to the wind. 
The slab of column density N env is at a single temperature T mv 
and moves at velocity V e „„. The widths of the slab and wind ab- 
sorption lines are 5 and 10 km s -1 respectively as suggested by the 
values found in Fig. [4] 

We performed two series of fits. In the first series the wind 
parameters were bracketed around the best values found by fitting 
the ISAAC data. In the second series, all parameters were given the 
freedom to take any physically valid value. The best fit in the first 
series is plotted in red in Fig.[TT]and the best fit in the second series 
is plotted in green. The parameters for both fits, the \ 2 of the fits, 
and their asymptotic 95% Confidence Interval are summarized in 
Table [5] We only provide the asymptotic errors as the fits are far 
from optimal (xt > 100). Both fits, which are based on a constant 
mass-loss rate wind model, fail to reproduce the complex shape of 
the 12 CO absorption profile but the fit with a mass-loss rate an or- 
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Figure 10. Best fit by the wind model. The upper graph shows the best model synthetic spectrum shifted by +2.0, the middle graph is the observed ISAAC 
transmission spectrum shifted by +1.0, and the lower graph shows the residual to the best fit. The best fit does not completely account for absorptions of very 
low- and very high- J transitions. 



der of magnitude lower than the value found by fitting the ISAAC 
data seems to give a better match. A closer inspection reveals that 
both fits struggle in the region between -30 and -50 km s~ . The 
low mass-loss rate model underpredicts the amount of absorption 
whereas the high mass-loss rate model overpredicts the absorption 
feature. The high mass-loss rate fit over-predict the amount of ab- 
sorption while the low mass-loss rate fit under-predict the amount 
of material in this velocity range. The mass-loss rates provided by 
the two fits bracket the mass-loss rates. The temperature of the slab 
is between 51 and 83 K, warmer than the value derived from the 
Gaussian fitting analysis. As expected from the 13 CO ISAAC data, 
the 12 CO column density is three orders of magnitude higher, be- 
cause the rotational diagram analysis vastly underestimates column 
densities in very optically thick cases. The temperature profile is 
not well constrained. Both a slightly increase or decreasing tem- 
perature profile are consistent with the data. 



7 DISCUSSION 

7.1 Gaseous and solid CO around LLN 19? 



From our analysis of the 
7V icc (CO)W gaa (CO) ratio is - 



gas and solid phase CO, the 
5 x 10~ 3 . The CO ice/gas is the 



lowest seen towards other YSOs in the Vela molecular cloud with 
similar bolometric luminosity (£, = 1600 L© for LLN 19) but close 
to the value derived in a few high-mass YSOs (L JLq >10 4 ) 
which have less than 1% of CO in the solid state dMitchell et al.l 
1990). Other high-mass YSOs have much higher fraction of CO 
in the solid state (i.e., 22.7 % CO ice in the envelope of RAFGL 
7009S jDartois et al.ll 19980 ). Also, contrary to other young stellar 
objects in the Vela cloud, LLN 19 has a low water ice abundance 
with respect to H2 of 1.1 x 10" J and a tot al water colu mn density 
of A icc (H 2 0)= 1.05 ± 0.2 x 10 18 cm -2 dThi et alfeOOfj) . One ex- 
planation is that LLN 19 is in a m ore evolved stage than the other 
YSOs studied in iThi et alj (2006) and has already cleared its sur- 
rounding envelope via powerful outflows. In this case, both wa- 
ter and CO ices are most likely located in the foreground material 
along the line-of-sight toward LLN 19. 

Another way to constrain the amount of gas phase CO deple- 
tion is to compare the estimated total column density CO molecules 
in the line-of-sight from extinction studies and the detected amount 
of CO in the gas phase. The total CO column density in the slab 
model is 8.9 x 10 18 cm" 2 and 3.5 x 10 18 cm" 2 in the wind model. 
In comparison, the fit to the 10 /J,m silicate feature usin g a ID ra- 
diativ e transfer code gives A(CO) ~ 8.5 x 10 18 cm -2 dThi et al.l 
120061) assuming [H 2 ]/[CO]=10 4 , i.e. without CO depletion. The 
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Figure 11. Two fits to the CRIRES observations. The red fit results from a wind+envelope model with mass loss rate of ~ 4.2 X 10~ 7 yr _1 while the green 
fit results from a model with mass loss rate of ~ 4.8 X 10 -8 yr -1 . The red fit grossly overestimates the absorption around the -43.5 km s" 1 features and fails 
to match the data in the blue tail at ~-70 km s _1 . On the other hand the green fit under- estimates the amount of absorption. Both fits stress the failure of the 
constant mass-loss rate model. 



Table 5. Best fit parameters to the CRIRES data for the two series. C. I. 
means confidence interval. The confidence intervals are approximate statis- 
tical errors (see Appendix). Systematic errors are not included. 



Units Best fit 95% C. I. Best fit 95% C. I. 

low m high m 





AU 


0.14 


± 0.02 


0.9 


±2.6 


Rout 


AU 


73 


±4 


167 


-167/+298 


T 


K 


290 


±2 


340 


±6 


p 




-0.19 


±0.08 


0.07 


±0.8 




km s — 1 


16.0 


±0.02 


15.9 


0.7 


Vmax 


km s — 1 


70 


±0.6 


28 


± 1 


<? 




2.2 


±0.08 


2.3 


± 1.2 


m 


10~ 7 M yr" 1 


0.48 


±0.02 


4.2 


±3.5 




10 17 cm" 2 


4.3 


±0.015 


6.0 


± 1.8 


Tslab 


K 


51 


± 1 


83 


±3 


Vslab 


km s — 1 


-8.9 


±0.03 


-8.6 


±0.02 






125 




224 





similarity between the observed and estimated CO column den- 
sity suggests that CO is lightly frozen onto grain surfaces along 
the line-of-sight towards LLN 19. 

We can also estimate the total CO column density from the 
visual extinction. We adopt the conversion factor Nh(CO)Mv ~ 
9 x 1 16 cm -2 mag -1 if there is no CO depletion (Frerki ng et al.l 
1982). The visual extinction is estimated from the observed E(H — 
K) with A v = (15.3 ± 0.6) x (H - K) jRieke & Lebofskvlll985h . 
The intrinsic H — K color of early B up to F stars is close to zero. 
Assuming that the K band flux is dominated by the extinct central 
source, E(H — K) ~H — K . The H- and if-ba nd magnitudes are 
11.92 and 8.88 respectively jLiseau et al. l [l992T) . which translates 
into an extinction of ~46 magnitude. The derived value of 7V(CO) 
is 4.2 x 10 18 cm -2 , smaller but close to the detected total value of 
8.9 x 10 ls cm -2 in the 2-slab model and very close to the value 
of 3.5 x 10 18 cm~ 2 found using the wind model. The observed 
gas-phase CO column density accounts for the predicted column, 
indicating that most CO molecules are in the gas phase. 

The two different fitting methods agree on the fact that LLN 19 
shows a much larger amount of warm than cold gas. The tempera- 
ture of the cooler component (45 K in the 2-slab model, 200-300 K 
in the wind model) is much higher than the sublimation temper- 
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ature of the pure CO ice (20 K) and is consistent with the low 
7Vicc(CO)/iVga a (CO) ratio of ~ 5 x 1(T 3 . Because LLN 19 is lo- 
cated in a dense stellar cluster environment, its outer envelope may 
be heated above the CO condensation/sublimation temperature by 
the ambient radiation field generated by the nearby stars. This phe- 
nomenon has been shown to work to explain molecular depletion in 
enve lopes around interm ediate-mass YSO dj0rgensen et al.| [2006). 

iRettig et alj (2006) found higher 7V ga s(CO) than expected 
from extinction studies of four YSOs and suggest that dust grains 
are settling in the disc midplane. The presence of an outflow may be 
considered an indirect evidence for the existence of a circumstellar 
disc around LLN 19. The high velocity attained by the outflowing 
gas and the absence of redshifted absorption indicate that the cir- 
cumstellar disc would be seen at low inclination with respect to the 
axis of rotation and the line-of-sight towards LLN 19 would inter- 
cept a small amount of disc material only. However, the extinction 
in the wind itself (Ay >50) is high enough to mask the redshifted 
absorption by recessing outflowing gas. 

A low CO depletion found in o ur study concu r s with obser- 
vations of the CS millimeter line by iGiannini et alj 120051 which 
suggest that no molecular depletion occurs in the Vela Molecular 
Cloud on a global scale. 



7.2 12 CO/ 13 CO ratio 



iGoto et alj d2003h compared optically thin 12 CO v = 2 <- 
and 13 CO v = 1 <— towards three YSOs. They derived a 
12 CO/ 13 CO ratio of 137 ± 9, 86 ± 49 and 1 58 toward LkHq 101 
AFGL 490, and Mon R2 IRS 3 respectively. IScoville et ail dl983h 
found a 12 CO/ 13 CO ratio of ~ 100 toward the BN/KL object. The 
ratios are 1.5-2 times higher than the generally accepted interstellar 
value. We estimated a 12 CO/ 13 CO ratio of 67 ± 3 towards LLN 19, 
close to the average interstellar value. The line-of-sights in Goto 
et al. study probe quiescent (dv=2-3.5 km s _1 ) cold foreground 
clouds whereas our data suggest that most of the CO gas is located 
in a turbulent warm region around LLN 19. 

More recently, a detail study of CO isotopologues absorption 
in the Af-band with CRIRES toward the protoplanetary disc VV 
CrA and the YSO Reipurfh 50 gives a value of C( 12 C/ 13 C) of 100 
dSmith et alj|2009h . One possible explanation for the higher ratio 
is selective photodissociation of 13 CO in protoplanetary discs and 
protostellar environment. It appears that the 12 CO/ 13 CO may be in- 
homogeneous within a cloud. This supports the idea that local dif- 
ferentiation processes occur. Detailed analysis of the 12 CO/ 13 CO 
ratio using optically thin v — 2 <— 12 CO absorption of the cold 
component at 9 km s _1 will be published elsewhere (Smith R. et 
al.). In terestingly, the C( 12 C/ 13 C) for CO ice has been found to be 
~ 70 jBoogert et al]|2002l ; Fontoppidan et alj|2003h , close to the 
quiescent interstellar gas phase values (60-80). 



7.3 Properties of the outflowing gas 

We derived some properties of the outflowing gas in LLN 19 and 
summarised them in Table[6] 

The geometry of absorbing gases cannot be constrained with 
spatially unresolved observations. The gas may be outflowing in 
form of spherical shells or in highly collimated jets. For simplicity, 
we have assumed that the wind is spherical, purely molecular, and 
that the n(CO)/n(H 2 ) ratio is 10" 4 . 

LLN 19 has a luminosity of 1600 L©. Its inferred mass-loss 
rate of 4.8 x 10~ 8 - 4.2 x 10~ 7 M yr -1 is two to three orders of 



magnitude lower than that in mid- to early-B (L« > 10 4 L ) YSOs 
~ 10 ~ 5 to a few 10" 3 M Q yr -1 dMitchell et al.ll 199ll : lArce et al.1 
l2007h . Low mass-loss rates of the order of 10 7 Mp yr 1 are f ound 
for low-mass (i* ^ 1 Lq) young stars dBontemps et al.l 1996b . Our 
result is margina lly consistent with a scaling of the mass-loss rate 
with luminosity dNisini et aljl995h . 

The source of the momentum carried by extended bipolar 
outflows remains undetermined. By computing the momentum of 
the young outflows we may compare them to values derived from 
millimeter observations of extended outflows. The assumption of 
spherical flow in our modelling implies that the derived momenta 
are lower limits. The momentum per unit area of a shell, assuming 
a thin shell compared to its radius R, is 



Psheli(-R) = m 



dR 



AtyR 2 ' 

and the momentum at the outer radius R ou t is 

(•Rout 

P — I mdR ~ mR out . 



(10) 



(11) 



The outer radius of the wind from the model is Rout — 73- 
480 AU (or (1-7.2) x 10 15 cm), which is consistent with values 
found for high-mass stars. The momentum per unit area in the wind 
at the outer radius is given by dMitchell etaljl99lh : 



p(Rout) = N{CO) 



n(H 2 ) 
n{CO) 



i{H 2 )V w = 



^TTRou 



(12) 



We derived a momentum per unit area p(R ut) of (2.9- 
8.4) x 10 2 g cirT 1 s _1 and a total momentum P of (0.19- 
9.6) x 10~ 4 M© km s" 1 . The momentum per unit area is ~ 100 
times and the total momentum is also ~ 1000 times smalle r than 
in outflows around high-mass YSOs dMitchell et a I] |l99lh . The 
momentum rate (also called momentum flux or force) of (0.4- 
7.9) x 10 -5 Mq km s _1 yr -1 is als o similar to the r ates found 
for outflows from low-mass protostars dArce et al]|2007h . 

We also estimated the dynamical timescale of the wind, which 
corresponds to the time for a gas shell to reach the wind outer radius 
Rout from Ri n . The mass in a shell at distance R from the source 
is: 



dm = 4nR mH 2 dN 



The constant mass-loss rate can be rewritten as 

. dm a dN(R) . p2 lV(flo„ 

m = — = AttR m H2 — — ^ 47r R out m H2 — 

at at tdyn 

The dynamical timescale is thus 

. 2 iV(fl out ) 

Idyn — 4TTR out m H2 : , 



(13) 



(14) 



(15) 



with TV = 10 4 x iV( 12 CO) = 3.5 x 10 22 cm" 2 at R out . We 
find that td, yn ~ 28 yr at the outer radius. Using the fits to the 
CRIRES data, we estimate the dynamical timescale to range be- 
tween 5 and 28 years. Therefore, whatever the actual value, the 
dynamical timescale is much shorter than the values derived from 
extended molecular outflows in the millimeter. 

The dynamic timescale is consistent with the values found for 
high-mass YSOs (2-200 yr). The total momentum and the dynam- 
ical timescale for the extended outflow obtained by mapping the 
environment of LLN 19 (8 ' x 4 ') i n 12 CO J = 1 — » is 17 1 
M km s" 1 and t dyn ~ 1.6 x 10 5 yr, IWouterloot & Brandl 19991) . 
Those values were derived assuming a distance of 2.24 kpc. The 
discrepancy between the millimeter and infrared values suggests 
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Figure 12. Schematics showing different outflow models. The outflow mo- 
mentum, temperature, and velocity profile along the axis are shown for the 
four major types of models: turbulent jet, jet bo w shock, wide-ang le wind, 
and circulation. The schematics are adapted from Arce et alj fe007l) . 



that the CO fundamental lines are probing the recent history (< few 
100 years) of the outflow. Outflows change morphology with time 
dArce et alj |20071. The youngest outflows tend to be highly colli- 
mated or include a very collimated component while older sources 
are less collimated with wider opening angles. Finally, the sharp 
drop in total momentum indicates that the outflow was probably 
much more active in the past. 

The high-resolution CRIRES data confirm the picture drawn 
from the fitting of the medium-resolution ISAAC data. The sub- 
structures seen in the CO absorption lines at different velocities fur- 
ther suggest that the outflow mass-loss is episodic with at least two 
major events. We bracketed the mass-loss rate between 0.48 and 4.2 
x 10 -7 Mq yr _1 . This is an order of magnitude variation in mass- 
loss and can be achieve d in the context of FU Orionis type burst 
fe.g. JCalvet et al Multiple blue-shifted absorptions are also 

found i n high-mass YSOs a nd are attributed to episodic mass-loss 
events i lMitchell et~aill99lb . 

Figure [T2l summarizes schematically the four major outflow 
models. From all the four models, only the jet bow shock model 
predicts episodic variation in jet properties as observed towards 
LLN 19 ( I Arce et all 1 2007b . Also the jet bow shock model is the 
sole model that can explain a temperature and velocity increase 
with distance from the outflow source. 

Another peculiar feature of the CO lines towards LLN 19 is 
the intrinsic line widths. Generally, 12 CO ro-vibrational absorp- 
tions/emissions toward YSOs are dominated by a hot component 
at T gaa > 600 K (e.g.. lMitchelletalJll990b . The intrinsic width 
is low 3-5 km s _1 in contrast to our result (d«=10-12 km s _1 ) 
for the outflowing gas, although widths up to 8-10 km s have 
been observed. Narrow cold gas found in the line-of-sight of high- 
mass YSOs is attributed to gas in the extended envelope. The wind 
model-envelope fit to the CRIRES data also suggests that a cold 
envelope may be present in the line-of-sight. 

One possibility for the warm gas to show large intrinsic width 
is that turbulence is at play. The detection of turbulent gas would be 



tantalising. The decay of turbulent eddies could provide the neces- 
sary heating to keep large amount of gas at ~ 100- 1000 K as the ra- 
diative heating from the central star decreases with distances. One 
possible origin for the turbulence are Kelvin-Helmholtz instabili- 
ties caus ed by the interaction of atomic jet entraining molecular 
material iDelam ater et alj2003b . 

Conversely the jet bow shock model of outflows provides a 
scenario to generate turbulent gases (see fig. I12l > . The atomic 
jets from the inner disc surfaces impact the ambient gas creating 
shocks. In turn, the high pressure gas between the shocks is ejected 
sideways. This high pressure gas interacts with unperturbed molec- 
ular gas, generating turbulences. Shocks will heat the gas to several 
hundred Kelvin, consistent with the gas temperatures of the differ- 
ent substructures. The jet bow model can explain the slight increase 
of gas temperature with radius contrary to the turbulent jet model. 
However the 95% confidence interval for the gas temperature al- 
lows decreasing gas temperature with distance to the star (Table[4]l. 
A flat temperature profile would be consistent with the wide-angle 
wind model. The large uncertainties in the model parameters, espe- 
cially for the gas temperature profile, prevent us to favour a specific 
outflow model. 



8 CONCLUSIONS 

We analysed medium- (7? ~ 10,000) and high-resolution 
(R ~ 100,000) 4.5—4.8 /im spectra of the embedded intermediate- 
mass young stellar object LLN 19 in the Vela Molecular Cloud. Gas 
phase 12 CO and ls CO ro-vibrational lines have been analysed with 
a curve-of-growth and a rotational diagram, and also fitted by a 
wind model. Both analyses give similar values for the total column 
density of warm CO gas and a 12 C/ 13 C isotopic ratio of 67 ± 3, 
lower than values close to 100 found in other studies using infrared 
absorption lines. The discrepancy in 12 C/ 13 C isotopic ratios de- 
rived from CO may just reflect the occurrence of varied chemical 
phenomena in molecular clouds and star-forming regions. 

The CO gas towards LLN 19 arises from an outflowing warm 
gas. Most CO molecules are in the gas phase with less than 0.5% 
condensed on grain surfaces, maybe in a foreground cloud. The 
wind parameters of LLN 19 are closer to those from low-mass pro- 
tostars than high-mass protostars, as seen in infrared absorption. 
The warm outflowing gas has been created only recently with a dy- 
namical timescale of less than 28 yr, if the data are analysed in the 
context of a static spherical wide-angle wind model. However, the 
kinematics of the outflow as seen in the high spectral resolution 
CRIRES data are more consistent with a episodic jet bow shock 
model where the internal bow shocks lead to high pressure and high 
temperature shells. A dynamical mass-loss rate model is warranted. 

LLN 19 (IRAS 08470-4321) is a good example of an object 
in the phase of cleaning its envelope by an outflow. Temporal mon- 
itoring of LLN 19 with CRIRES would provide further clues on the 
evolution of its outflow. 
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Table 6. Outflow parameters. 





Units 


ISAAC 


CRIRES 
low fa 


CRIRES 
high fa 


Total column density iV(H2) 


10 22 cm" 2 


3.5 


2.2 


2.6 


Mass loss rate fa 


1CT 7 M yr" 1 


4.2 


0.48 


4.2 


Momentum per unit area p 


10 2 g cm -1 s _1 


2.9 


2.5 


8.4 


Total momentum P 


M© km s _1 


9.6 x ur 4 


1.9 x ur 5 


3.3 x 1(T 4 


Momentum rate 


M Q km s _1 yr -1 


7.9 x 10" 5 


3.9 x 10" 6 


1.2 x 10- 5 


Dynamical timescale t dyn 


yr 


26 


5 


28 
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10 APPENDIX 

For nonlinear models there is no exact way to obtain the covariance 
matrix for the parameters. Nor do the optimization codes give di- 
rectly error parameters. Parameter confidence intervals are approx- 
imated from the Hessian matrix of T rno d(9) at the s olution 9 = 9 
jPonaldson &~Sc hnabel 1987; Wall & Jenkins 200l): 



H(9) = 2J(9) T J(9) + 2 fi (I) , 



(16) 



3=1 



where J (9) is the Jacobian of T mo< i(p), In the neighbourhood of the 
solution, \\fj(9)Hj(9)\\ is neglectable compared to \\J(9) T J{9)\\. 
Therefore, the Hessian matrix at the minimum for F{9) (equal to 
the residual sum square RSS) at 6* = 9 can be approximated by 

H{9) ~ 2J{9) T J{9) (17) 

An unbiased estimate of the parameter covariance matrix is 



C(9) 



(18) 



where s 2 = F(9)/dof is an unbiased approximation of the resid- 
ual variance a (a ~ s). The variance of the ith parameter var(9i) 
is the ith diagonal element of the matrix C(9). If 9* is the "true" 
solution, then the 100(1 — <f>) percentage confidence interval on 9t 



(19) 



where t 



l-<t>/2 
dof 



is the 100(1 — 4>/2) percentage point of the two- 
tail student t-distribution with dof degrees of freedom. When the 
number of degrees of freedom exceeds 150, the Student distribu- 
tion equals the normal distribution. In our case, a 95% confidence 
interval would be the parameter ±1.96 x \J var0i). The Jacobian 
matrix element Jij (9) is approximated by forward finite difference 
around 9 



Tmod^j, 9i) — T mo d(^j, 9i) 



(20) 



In this paper, the parameter intervals are computed for 95% confi- 
dence. 
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